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(57) Abstract: A new and improved method for tracking and/or 
spatial localization of an invasive device in Magnetic Resonance 
Imaging (MRI) is provided. The invention includes providing an 
invasive device including a marker having a chemically shifted 
signal source with a resonant frequency different from the chem- 
ical species of the subject to be imaged, applying a pulse se- 
quence, detecting the resulting RF magnetic resonance signals, 
and determining the 3D coordinates of the marker. The inven- 
tion also includes generating scan planes and reconstructing an 
image from the detected signals to generate an image having 
the marker contrasted from the subject. The invasive device in- 
cludes a marker having a chemically shifted signal source which 
has a resonant frequency different from the chemical species of 
the subject to be imaged for use in tracking the device during 
imaging. 
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CHEMICAL SHIFT MARKERS FOR IMPROVED 
WIRELESS FIDUCIAL MARKER TRACKING 

BACKGROUND OF THE INVENTION 

The present invention relates to a system and method for spatial localization of 
an invasive device in Magnetic Resonance Imaging (MRI). More particularly, the 
5 invention relates to a system and method for tracking an invasive device during MRI 
imaging. 

MRI provides visual information about the interior of a subject which can be 
very useful for surgical procedures. Using MRI, surgical tools such as interventional 
devices, also known as invasive devices, can be used to minimize the size of incisions 

10 resulting in shorter recovery times. The use of small incisions requires the invasive 
devices to be guided or tracked since they cannot be directly observed within the 
subject. MR-guided device tracking has generated considerable interest as a result. 

One key component of current tracking methods is the ability to automatically 
define an appropriate MRI imaging scan plane. Optical navigational systems have 

15 been developed to provide interactive scan plane selection by measuring the position 
of reflective materials mounted onto the surface of a rigid interventional device. 
Optical tracking systems have limitations including requiring a line-of-sight between 
the markers and the detection system. This hinders positioning of the detection 
systems within the scan room which in turn restricts the range-of-motion of the 

20 physician during the procedure. Also, these optical detection systems add additional, 
and potentially expensive, components to an already congested physical environment 
and they require calibration between the tracking system and the MRI coordinate 
frames. 

An improved wireless tracking method was recently developed using 
25 inductively coupled tuned fiducial markers, a limited projection reconstruction 
sequence (LPR-FISP), and a fast localization algorithm to provide automatic scan 
plane selection for interventional procedures. This tracking method improved upon 
optical tracking systems by providing fast, automatic scan plane selection without a 
line-of-sight requirement between the markers and the detection system. 
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Unfortunately, this system of fiducial marker contrast using tuned coils is dependent 
on the position and orientation of the marker coil with respect to the transmitter and 
receiver coils. In certain orientations, this dependence can result in attenuation of the 
marker signal and a loss of tracking function. It is desirable to provide a wireless 
5 tracking system and method to track markers which reduces the dependency of the 
signal strength on the position and/or orientation of the marker. 

SUMMARYOF THE INVENTION 
According to the present invention, a new and improved method for spatial 
1 0 localization of an invasive device in Magnetic Resonance Imaging (MRI) is provided. 

In accordance with a first aspect of the invention, the invention includes 
providing an invasive device including a marker having a chemically shifted signal 
source with a resonant frequency different from the chemical species of the subject to 
be imaged, applying a pulse sequence, detecting the resulting RF magnetic resonance 
1 5 signals, determining the 3D coordinates of the marker. 

In accordance with second aspect of the invention, the invention includes 
selecting scan planes for MR imaging. 

In accordance with another aspect of the invention, the invention includes 
reconstructing an image from the detected signals to generate an image having the 
20 marker contrasted from the subj ect. 

In accordance with a yet another aspect of the invention, an invasive device 
for use in Magnetic Resonance Imaging (MRI) of a subject is disclosed having a 
marker with a chemically shifted signal source which has a resonant frequency 
different from the chemical species of the subject to be imaged for use in tracking the 
25 device during imaging. 

In accordance with a yet another aspect of the invention, a method for 
selecting a scan plane in MRI is provided. The method includes providing an 
invasive device including a marker having a chemically shifted signal source with a 
resonant frequency different from the chemical species of a subject to be imaged, 
30 applying a pulse sequence, detecting the resulting RF magnetic resonance signals, 
determining the 3D coordinates of the marker, and selecting a scan plane which 
include the 3D coordinates of the marker. 
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Other features, benefits and advantages of this invention will become apparent 
to those skilled in the art from the following detailed description of the preferred 
embodiments, when read in light of the accompanying drawings. 

5 

BRIEF DESCRIPTION OF THE DRAWINGS 
The invention may take form in certain components and structures, preferred 
embodiments of which will be illustrated in the accompanying drawings wherein: 

10 Fig. 1 illustrates a method of invasive device tracking in accordance with the 

invention; 

Fig. 2 is an illustration of an invasive device with markers attached thereto in 
accordance with the invention; 

Fig. 3 is a radial 2-pulse missing pulse steady state sequence; 
1 5 Fig. 4 illustrates steps of the invention; 

Fig. 5 illustrates projection signal peaks used in determining the 3D 
coordinates of the fiducial markers in accordance with the invention; 

Fig. 6 illustrates steps of the invention; 

Fig. 7 illustrates the intersections created by backprojecting the marker signal 
20 peaks from each projection in accordance with the invention; 
Fig. 8 illustrates steps of the invention; 

Fig. 9 illustrates closest-point sets used for determining the 3D coordinates of 
the fiducial markers in accordance with the invention; 

Fig. 10 illustrates steps of the invention; 
25 Fig. 11 illustrates centroids, of the N= 3 densest closest-point sets used to 

represent the actual 2D locations of the fiducial markers in the scan plane; 

Fig. 12 illustrates steps of the invention; 

Fig. 13 illustrates the matching of the corresponding Y coordinates of the 
fiducial markers from axial (X-Y plane) and sagittal (Y-Z plane) scans in accordance 
30 with the invention; and 

Fig. 14 illustrates a method for selecting a scan plane in accordance with the 
invention. 
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DETAILED DESCRIPTION OF THE INVENTION 
Referring now to Fig. 1, a method for spatial localization of an invasive device 
in a subject using MRI is shown generally at 10. The method 10 includes providing 
an invasive device at 12 with a marker, also known as a fiducial marker, having a 
5 chemically shifted signal source which has a resonant frequency different from the 
chemical species of the subject to be imaged. The fiducial markers provide a set of 
highly localized signal sources that are tracked to determine the location and/or 
orientation of the invasive device. 

The invasive device can have one or a plurality of markers attached thereto, 
10 disposed either in, on, or about the device. Referring now to Fig. 2, an invasive 
device D is shown having three markers 13 attached to the outer surface of the device. 
The invasive device D can be a catheter, a guide wire, an endoscope, a laparoscope, a 
biopsy needle, a surgical instrument or any other suitable known interventional or 
invasive device. 

15 The markers 13 each have a chemically shifted signal source with a resonant 

frequency different from the chemical species of the subject to be imaged (not 
shown). Markers were created by infusing 2.5ml of concentrated acetic acid (0 ~ 
lOppm)) doped with ImM Gd contrast (OPTMARK™, Mailinckrodt Inc.) into 12 ml 
syringes (ID = 15mm) to provide a signal source with a proton chemical shift distinct 

20 from tissue protons encountered in a typical scanned volume. However, any suitable 
marker having any suitable amount of a chemically shifted signal source different 
from the chemical species of the subject to be imaged can be used. Other examples, 
which should not be considered limiting, can include blood-safe compounds including 
fluorinated compounds, such as perflourocarbon or any suitable known materials 

25 having a unique resonant frequency different from the chemical species of the subject 
to be imaged. 

The one or more markers 13 can be disposed in any suitable known 
arrangement. The use of a known arrangement enables the one or more markers 13 to 
provide device location and/or orientation information. For example, the marker 13 
30 can be arranged as a point source. A single marker 13 can be disposed at the tip of the 
device D to indicate the location of the device tip. The marker 13 can be arranged as 
a line source, with the line providing device orientation information. Several markers 
13 can be used as points sources. The point sources can be arranged in a line, or three 
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or more can be arranged to form a plane as shown in Fig. 2. The plane can provide 
useful device orientation information. 

Referring again to Fig. 1, the method 10 further includes applying an MR 
pulse sequence to at least a portion of the subject at 14. A fast, radial MRI pulse 
5 sequence was used for detecting the local signal from the fiducial markers and 
suppressing other signals in the image volume. 

Referring to Fig. 3, a radial 2-pulse missing pulse steady state (MPSSFP) 
sequence (TR/TE/FA=30ms/20ms/135°) shown generally at 30 was used. The pulse 
sequence includes a slice-selective RF excitation pulse (TH=100mm) shown at 32 and 
10 a chemical shift selective RF excitation pulse shown at 34 to generate marker-only 
projections. If scan volumes are used, the pulse 32 is a volume selective RF 
excitation pulse. 

Dephasing as shown at 36 was applied in the readout direction following the 
slice-selective pulse 32 to remove residual signal from the subject. The slice-selective 

15 RF excitation 32 was used to prevent off-resonance spins at the edge of the main 
magnetic field from contributing to the marker-only images. Alternatively, any 
known pulse sequence suitable for device tracking can be used. 

The method 10 further includes detecting the resulting RF magnetic resonance 
signals at 16. The signals are detected in any suitable manner. Scan planes 

20 containing the marker or markers 13 are acquired. The pulse sequence can generate a 
limited set of ID projections that can be analyzed by the localization algorithm as 
described below. Ideally, the sequence would detect only projection information from 
the fiducial markers while suppressing signal from all other sources (e.g., patient or 
water phantom, here). Alternatively, scan volumes containing the marker or markers 

25 13 can be acquired in any suitable known maimer. 

The method 10 further includes determining the 3D coordinates of the marker 
or markers at 1 8. A computer algorithm was used to analyze the ID projections from 
the pulse sequence and to identify the 3D coordinates of the fiducial markers 13. 

A five-stage method to accurately determine the 3D coordinates of the fiducial 

30 markers 13 was developed. However, this method is given for the purposes of 
example only and the 3D coordinates of the fiducial markers can be determined in any 
suitable known manner from the scan planes or scan volumes. In phase one, the 
algorithm converts the raw MRI projection data into the spatial domain, by applying a 
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1D-IFFT, and finds the location of the markers within each projection. In phase two, 
the algorithm analytically determines the location of all the intersections that would 
be created by backprojecting the marker signal peaks that were found in each 
projection into the scan plane. Phase three involves designating a subset of these 
5 intersection points as reference points and forming a "closest-point set" around each 
of these points. In phase four, the centroids of the N densest "closest-point sets" are 
used to represent the 2D locations of the fiducial markers in the scan plane (with N 
equal to the number of markers used, in this example N = 3). Finally, in phase five, 
the 3D locations of the markers are calculated. The XY and YZ coordinates are found 

10 separately using the aforementioned method, and by matching the corresponding Y 
coordinates, 3D locations are assigned to each fiducial marker. 

Referring now to Fig. 4, the step of determining the 3D coordinates of fiducial 
markers 18 includes converting the raw MRI projection data into the spatial domain, 
by applying a 1D-IFFT, and finding the location of the markers within each projection 

15 shown generally at 40. The algorithm was designed to perform marker localization 
using a pre-determined number of projections in each scan plane (typically this value 
was set to 5, although any suitable number may be used). For a projection to be 
useful to the localization algorithm, all of markers signals should be identified within 
the projection. In this example, 3 markers were used, although any suitable number 

20 of markers can be used. 

The raw MRI projection data was converted into the spatial domain by 
applying a 1D-IFFT at 42. Numerical algorithms were then used to search for fiducial 
marker signals in the projections by identifying the 3 largest signal peaks at 44. A 
discrete differentiation formula was used to identify the 3 largest local maxima in 

25 each projection. For the purposes of example, the 3 peaks in a projection are shown at 
45 in Fig 5. 

Next it was verified that all of the peaks 45 met an experimentally determined 
minimum signal to noise ratio (SNR) threshold at 46. The SNR of each peak was 
then calculated by dividing the intensity of the signal peak by the mean background 
30 noise level in the projection. This noise level was estimated by removing all 3 marker 
signal candidates (which included the signal peaks and four data points on either side 
of each peak) from the projection data. If any of the 3 marker signal candidates did 
not meet the minimum SNR requirement, the entire projection was discarded at 48. 
This step was continued until enough satisfactory projections were identified at 50. 
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Referring now to Fig. 6, the step of determining the 3D coordinates of fiducial 
markers 13 also includes analytically calculating projection intersections as shown 
generally at 60. The algorithm determines at 62 the location of all the intersections 63 
that would be created by backprojecting the marker signal peaks 45 from each 
5 projection as shown in Fig. 7. 

This intersections were determined analytically by representing every marker 
signal peak 45 by a linear equation defined by a slope (tan9) and an intercept (S/cos0). 
Here 0 is the projection angle and S is the location of the peak along the projection 
axis. The solution of any two linear equations (provided they do not correspond to 
10 marker signals from the same projection) defines an intersection point 63 in the scan 
plane. The linear equations from a pair of projections define N 2 intersection points, 
and N of these points corresponds to the location of a fiducial marker (with N = 3 in 
this example). By solving for the intersections generated by every possible 

r-i 

combination of projections, L 2 J *N 2 intersections were found in the scan plane. Here 
15 m is the number of projections used by the algorithm. The projections were ranked 
according to their d m i D as described below. 

Referring now to Fig. 8, the step of determining the 3D coordinates of fiducial 
markers also includes generating Reference Points and Closest-Point Sets shown 
generally at 70. Intersection points generated by some pair of projections were 
20 designated as reference points in the scan plane at 72. The two projections used to 
define the reference points were those that, had the largest distance separating their 
two closest marker signal peaks 45, d m | n as shown in Fig. 7, to ensure that the 
reference points were well distributed in the scan plane and not too close together. A 
"closest-point set", shown as 76 in Fig. 9, containing the intersection point 63 from 

25 each of the L 2 J projection pairs that was closest to the reference point, was formed 
around each reference point at 74. Closest-point sets were generated by iteratively 
searching through each set of intersection points and selecting the point that was 
geometrically closest to a given reference point. 

Referring now to Fig. 10, the step of determining the 3D coordinates of 

30 fiducial markers also includes estimating the closest-point set density and centroid. 
The mean separating distance was computed for each closest-point set by averaging 
the distances between the reference point and the other members of the set. A closest- 
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point set's relative density was estimated by comparing the mean separating distances. 
Each set point was ranked according to its mean separating distance at 82. The 
centroids, shown at 83 in Fig. 11, of the N= 3 densest (smallest mean separating 
distance) closest-point sets were used to represent the actual 2D locations of the 
5 fiducial markers in the scan plane at 84. 

Referring now to Fig. 12, the step of determining the 3D coordinates of 
fiducial markers also includes obtaining the 3D coordinates using orthogonal scan 
planes shown generally at 90. A linear least squares routine was used to match the 
corresponding Y coordinates of the fiducial markers from axial (X-Y plane) and 
10 sagittal (Y-Z plane) scans as shown in Fig. 13. The resulting three-dimensional 
position of each fiducial marker 13 included the average of the two Y coordinates, 
and the X and Z coordinates that were associated with the marker's 2D position in 
each scan plane. 

The method 10 also includes reconstructing an image from the detected 

15 signals, shown at 20 in Fig. 1, in any suitable know manner to produce an image of 
the marker for use in tracking the device. The radial MPSSFP images were 
reconstructed from 256 projections, though any suitable number can be used. 

In accordance with a yet another aspect of the invention, a method for 
selecting scan planes in MRI is provided as shown generally at 100. The method 

20 includes providing an invasive device including a marker having a chemically shifted 
signal source with a resonant frequency different from the chemical species of a 
subject to be imaged at 102. The device can be any suitable invasive device as 
described above having one or more markers as described above. 

The method also includes applying a pulse sequence at 104 similar to the pulse 

25 sequence described above, and detecting the resulting RF magnetic resonance signals 
at 106. The method also includes determining the 3D coordinates of the marker at 
108 in any suitable manner, such as for example that described above. The method 
also includes selecting one or more scan planes which include the 3D coordinates of 
the marker or markers at 110. The method can also include reconstructing an image 

30 having the marker contrasted from the subject at 1 12. 

The invention improves tracking of markers 13 based on chemical shift which 
can be used to track and guide an invasive device D when attached thereto. The 
markers 13 can be disposed on, in or about the invasive device D to provide spatial 
localization of the markers which can be used for device tracking and acquisition of 
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scan plans. With the invention the position and orientation of the interventional 
device D is less constrained as compared to other known fiducial marker tracking 
methods. 

A radial missing pulse steady state free precession (MPSSFP) sequence can be 
5 used to provide marker contrast necessary for 3D localization of the markers and 
automatic scan plane selection in interventional MRI procedures, although other 
suitable known pulse sequences can be used. The marker contrast provided by the 
radial MPSSFP sequence allows the 3D coordinates to be analytically determined by 
existing localization algorithms. The radial MPSSFP was established to match the 
10 signal characteristics and relaxation parameters of the concentrated acetic acid / Gd- 
contrast solution. 

The invention has been described with reference to preferred embodiments. 
Obviously, modifications and alterations will occur to others upon reading and 
understanding the preceding specification. It is intended that the invention be 
15 construed as including all such modifications and alterations insofar as they come 
within the scope of the appended claims or the equivalents thereof. 
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1. A method for spatial localization of an invasive device in Magnetic 
Resonance Imaging (MBS) comprising: 

5 providing an invasive device including a marker having a chemically shifted 

signal source with a resonant frequency different from the chemical species of the 
subject to be imaged; 

applying a pulse sequence; 

detecting the resulting RF magnetic resonance signals; and 
1 0 determining the 3D coordinates of the marker. 

2. The method for spatial localization defined in claim 1 further 
comprising selecting scan planes for MR imaging. 

15 3 . The method for spatial localization defined in claim 2 wherein the scan 

planes include the marker. 

4. The method defined in claim 1 further comprising acquiring scan 
volumes containing the marker. 

20 

5. The method for spatial localization defined in claim 1 further 
comprising generating marker only projections for tracking the device during it's 
movement in the subj ect. 

25 6. The method for spatial localization defined in claim 1 wherein the 

marker is doped to improve the signal received in the detecting step. 

7. The method for spatial localization defined in claim 1 wherein the 
marker includes acetic acid. 

30 

8. The method for spatial localization defined in claim 1 wherein the 
marker is blood-safe. 
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9. The method for spatial localization defined in claim 8 wherein the 
marker includes fluorinated compounds. 



10. The method for spatial localization defined in claim 9 wherein the 
5 marker includes perflourocarbon. 

11. The method for spatial localization defined in claim 1 wherein the 
marker is arranged as a line source. 

10 12. The method for spatial localization defined in claim 1 wherein the 

device includes a plurality of the markers distributed as point sources. 

13. The method for spatial localization defined in claim 12 wherein the 
markers are arranged in a line to provide an image having a line defining the device 

1 5 contrasted from the subj ect . 

14. The method for spatial localization defined in claim 1 wherein the 
device includes at least three markers distributed as point sources defining a plane. 

20 15. The method for spatial localization defined in claim 14 further 

comprising using the plane to indicate the orientation of the device. 

16. The method for spatial localization defined in claim 1 wherein the step 
of applying a pulse sequence includes applying a radial 2-pulse missing pulse steady 

25 state sequence. 

17. The method for spatial localization defined in claim 16 wherein the 
radial 2-pulse missing pulse steady state sequence includes applying a slice-selective 
pulse and a chemical shift selective pulse. 

30 

18. The method for spatial localization defined in claim 16 wherein the 
radial 2-pulse missing pulse steady state sequence includes 
TR/TE/FE=30ms/20ms/l 35°. 
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19. The method for spatial localization defined in claim 1 wherein the step 
of applying a pulse sequence includes applying a dephasing pulse. 



20. The method for spatial localization defined in claim 1 further 
5 comprising reconstructing an image from the detected signals to acquire an image 

having the marker contrasted from the subject. 

21. The method defined in claim 1 further comprising generating a limited 
set of ID radial projections in two orthogonal scan planes. 

10 

22. The method defined in claim 21 wherein the determining step includes 
analyzing the ID projections from the pulse sequence to identify the 3D coordinates 
of the markers. 

15 23 . The method defined in claim 22 wherein the determining step includes : 

converting the raw MRI projection data into the spatial domain by applying a 
lD-FFt; and 

finding the location of the markers within each projection. 

20 24. The method defined in claim 23 wherein the determining step further 

includes locating the intersections created by backprojecting the marker signal peaks 
found in each projection into the scan plane. 

25. The method defined in claim 24 wherein the determining step further 
25 includes designating a subset of the intersection points as reference points and 

forming a "closest-point set" around the reference points. 

26. The method defined in claim 25 wherein the determining step further 
includes representing the 2D locations of the markers in the scan plane using the 

30 centroids of the densest "closest-point sets". 

27. The method defined in claim 26 wherein the determining step further 
includes finding the XY and YZ coordinates separately. 
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28, The method defined in claim 27 wherein the determining step further 
includes matching the corresponding Y coordinates to assign the 3D locations to each 
marker. 

5 29. The method defined in claim 1 further including tracking the invasive 

device. 

30. An invasive device for use in Magnetic Resonance Imaging (MRI) of a 
subject comprising: 

10 a marker having a chemically shifted signal source which has a resonant 

frequency different from the chemical species of the subject to be imaged for use in 
tracking the device during imaging. 

31. The invasive device defined in claim 30 wherein the marker is doped 
1 5 to improve the signal received during imaging. 

32. The invasive device defined in claim 30 wherein the marker includes 
acetic acid. 

20 33. The invasive device defined in claim 30 wherein the marker is blood- 

safe. 

34. The invasive device defined in claim 30 wherein the marker includes 
fluorinated compounds. 

25 

35. The invasive device defined in claim 34 wherein the marker includes 
perflourocarbon. 

36. The invasive device defined in claim 30 wherein the marker is 
30 arranged as a line source to provide an image having a line defining the device 

contrasted from the subject. 

37. The invasive device defined in claim 30 further including a plurality of 
the markers distributed as point sources about the device. 
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38. The invasive device defined in claim 30 wherein at least three markers 
are distributed as point sources defining a plane. 

5 39. The invasive device defined in claim 38 wherein the plane indicates 

the orientation of the device. 

40. A method for selecting a scan plane in Magnetic Resonance Imaging 
(MRI) comprising: 

10 providing an invasive device including a marker having a chemically shifted 

signal source with a resonant frequency different from the chemical species of a 
subject to be imaged; 

applying a pulse sequence; 

detecting the resulting RF magnetic resonance signals; 
1 5 determining the 3D coordinates of the marker; and 

selecting a scan plane which includes the 3D coordinates of the marker. 

41. The method defined in claim 40 wherein the invasive device includes a 
plurality of markers. 

20 

42. The method defined in claim 41 wherein the markers are arranged in a 

line. 

43. The method defined in claim 41 wherein the device includes at least 
25 three markers distributed as point sources defining a plane. 

44. The method defined in claim 43 further comprising using the plane to 
indicate the orientation of the device. 

30 45. The method defined in claim 40 wherein the step of applying a pulse 

sequence includes applying a radial 2-pulse missing pulse steady state sequence. 
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46. The method defined in claim 45 wherein the radial 2-pulse missing 
pulse steady state sequence includes applying a slice-selective pulse and a chemical 
shift selective pulse. 

5 47. The method defined in claim 46 wherein the radial 2-pulse missing 

pulse steady state sequence includes TR/TE/FE=30ms/20ms/135°. 

48. The method defined in claim 40 wherein the step of applying a pulse 
sequence includes applying a dephasing pulse. 

10 

49. The method defined in claim 40 further comprising reconstructing an 
image from the detected signals to acquire an image having the marker contrasted 
from the subject. 
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